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ABSTRACT 

We report on sensitive sub-mm imaging observations of the prototype 
Seyfert 2/starburst galaxy NGC 1068 at 850 /im and 450 /im using the 
Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk 
Maxwell Telescope (JCMT). We find clear evidence of dust emission associated 
with the extended HI component which together with the very faint 12 CO J=l-0 
emission give a gas-to-dust ratio of M gas /M dust ~ 70 — 150. This contrasts with 
the larger ratio M gas /M dust ~ 330 estimated within a galactocentric radius of 
r < 1.36 kpc, where the gas is mostly molecular and starburst activity occurs. 
The large gas-to-dust ratio found for the starburst region is attributed to a 
systematic overestimate of the molecular gas mass in starburst environments 
when the luminosity of the 12 CO J=l-0 line and a standard galactic conversion 
factor is used. On the other hand sub-mm imaging proves to be a more powerful 
tool than conventional CO imaging for revealing the properties of the diffuse H 2 
that coexists with HI. This molecular gas phase is characterized by low densities 
(n(H 2 ) < 10 3 cm -3 ), very faint emission from sub-thermally excited CO, and 
contains more mass than HI, namely M(H 2 )/M(HI) ~ 5. 

Subject headings: ISM: dust, molecules and atoms — galaxies: individual 
(NGC 1068) 
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1. Introduction 

The standard way to estimate molecular gas mass in galaxies uses the luminosity of 
the 12 CO J =1-0 line and converts it to molecular gas mass by using the so-called standard 
galactic conversion factor X co (e.g. [Young fc Scoville 1982[ [Bloemen 1985| ; |Dickman et al 



1986| ; [Young fc Scoville 1991|) . The dependence of X co on the ambient conditions of the 



molecular gas has been extensively explored. The main factors are: a) metalicity and the 
intensity of the ambient UV radiation field (e.g. [Israel 19"8"H| ; |Wilson 1995| ; |Arimoto, Sofue 



fc Tsujimoto 1996] ; [Israel 199"7| ), b) the density, temperature, and kinematic state of the 



average molecular cloud and, c) effects of particular geometries (e.g. see [Bryant & Scoville] 
1996[ |Sakamoto 1996| for a recent exposition). In the environments of extreme starbursts 



it has been clearly demonstrated that the standard method overestimates M(H 2 ) since it 



yields masses comparable or larger than the dynamical ones (e.g. Pownes, Solomon 



[Radford 1991$ |Bryant fc Scoville 1996t |Solomon et al. 1997| ; |Downes fc Solomon 1998|) . In 



low metallicity environments with high UV intensities the larger rate of CO dissociation 
with respect to H 2 produces a lower CO luminosity per H 2 and thus underestimates M(H 2 ). 

In principle the aforementioned effects can be taken into account by a combination of 
observational and theoretical work studying a wide range of galactic environments, which 
can then provide "adjusted" X co factors to be used for the appropriate set of conditions. 
However little can be done if the H 2 is so diffuse and/or cold that the 12 CO J=l-0 line is 
very faint or not not luminous at all. Hints of such physical conditions have been recently 
reported for the outer parts of Milky Way (Usuda et al. 1998) and cold dense clouds have 
been reported for M 31 ([Loinard fc Allen 1998|) . In such cases the advantage of easily 



detecting and mapping a bright CO line and thus the distribution of H 2 is lost. Here we 
report deep sub- mm imaging of the prototype Sy2/starburst galaxy NGC 1068 that reveals 
extended emission from dust well beyond the CO-bright regions and associated with the HI. 
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2. Observations 

The observations were made on two nights in 1997 August 4 and 1998 January 26 with 
the Sub-mm Common User Bolometer Array (SCUBA) at the 15 m James Clerk Maxwell 
Telescope (JCMT)Q. SCUBA is a dual camera system cooled well below 1 K that allows 
sky-background limited simultaneous observations with two arrays. The short- wavelength 
array contains 91 pixels and the long-wavelength array 37 pixels. Both arrays have 
approximately the same field of view on the sky ~ 2.3'. For a full description of the 
instrument see |Holland et al. 1998| . 

We performed dual wavelength imaging at 450 /xm and 850 /xm using the 64-point jiggle 
mapping mode that allows Nyquist sampling of the field of view ( [Holland et al. 1998| ). We 



employed the recommended rapid beam switching at a frequency of 8 Hz and a beam throw 
of 120" in azimuth. The pointing and focus of the telescope were monitored frequently 
using Mars, Jupiter and CRL 618, with an expected rms pointing error of < 3". All the 
NGC 1068 maps were bracketed by skydips used to correct for atmospheric extinction. 
Typical opacities at 850 /xm were r ~ 0.14 for our 1998 run and r ~ 0.6 for our 1997 run. 
Frequent photometric measurements and beam maps of CRL 618 and Mars allowed close 
monitoring of the gains (Jy/beam Volt -1 ). These can vary significantly especially at 450 /xm 
if the dish has not thermally relaxed and can even be elevation dependent (Sandell 1998). 

The individual jiggle maps are consistent in terms of peak and integrated intensities 
and were co-added after flatfielding, correcting for atmospheric extinction, and 
editing out bad bolometers/integrations using the standard reduction package SURF 



x The JCMT is operated by the Joint Astronomy Center in Hilo, Hawaii on behalf of 
the parent organizations PPARC in the United Kingdom, the National Research Council of 
Canada and the The Netherlands Organization for Scientific Research. 
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fljcnncss fc Lightfoot 1998 ). Special care was taken to remove sky-noise using bolometers 
at the edge of the field of view that "look" at sky emission only (|Jenness, Lightfoot fc 



[Holland 1998|) . As a result the final maps have exceptional sensitivity, close to the one 
expected from the total integration time and the NEFDs at the sky conditions of our runs. 
Flux calibration of the final maps in mJy/beam was obtained from beam maps of CRL 618 
obtained with the same chop throw as our NGC 1068 maps. The corresponding beam 
profiles were obtained from beam maps of CRL 618 and Mars. The calibration information 
is summarized in Table 1. 



EDITOR: PLACE TABLE 1 HERE 



3. Extended sub-mm emission: CO, HI and dust 

The clear association of the sub-mm emission from dust and HI gas is demonstrated 
in Figure 1 where the 850 fim and 450 /zm maps are overlayed with an HI map (Brinks 
et al. 1994) at a common resolution of ~ 15". Figure 2 shows the bright central 450 /j,m 
emission and the associated 12 CO J=l-0 integrated brightness ( [Heller fc Blitz 1995 ) that 



lie within the central HI "hole" seen in Figure 1. This wealth of imaging data for the ISM 
in NGC 1068 and the wide range of conditions present due to the existence of the central 
starburst (Telesco et al. 1984) make it an ideal testing "ground" of our standard ideas 
about the distribution and mass of the various ISM phases on large scales. We adopt a 
distance to NGC 1068 of 14 Mpc (H D = 75 km s" 1 Mpc" 1 ) where 1" = 68 pc. 

A remarkable property of the sub-mm emission is that it extends over the entire 
bright HI emission where the estimated R = S450/S850 ratio is lower than in the CO-bright 
region where the starburst occurs. This ratio and its variations across the image can be a 
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good temperature indicator as long as the average dust temperature is T < 30 K. Indeed, 
assuming a single average dust component, this ratio can be expressed as 



'850 



3 31.8/T _ I 



where (3 is the emissivity law. It can be easily seen that the ratio tends to the R-J limit 
for dust temperatures of T > 30K and is no longer temperature-sensitive. Nevertheless, 
for a range of T = 10 — 30 K, it varies by a factor of ~ 2. This temperature range is 
particularly interesting since it "marks" the transition between cool and warm dust, the 
latter dominating the IRAS 100 /xm and 60 /im bands. 

After carefully correcting the integrated flux densities for the error beam 
contribution (Sandell 1998) using the correction factors tabulated in Table 1, we 
obtain R(r < 20") = 7.10 ± 1.77. For (3 = 2 this agrees well with the range of gas 
temperatures T ~ 20 — 30 K found from studying the CO excitation within this region 
( Papadopoulos fc Seaquist 1998Q and the minimum temperature of ~ 20 K implied from 
the observed 12 CO J=l-0 brightness temperatures in high resolution maps (e.g. |Planesas~. 



Scoville, fc Myers 1991|) . For the sub- mm emission associated with the HI distribution 
we find R(20" < r < 60") = 4.70 ± 1.17. This ratio suggests colder gas, in the range 
of T ~ 10 — 15 K, but rules out any M 31-type clouds (Loinard & Allen 1998) where 
Tkin < 10 K. Here we must emphasize that the main uncertainty in the estimate of R is 
due to systematic factors (~ 25%) that stay unchanged across the maps. The thermal rms 
uncertainty is < 5% and thus the observed change of R between the CO-bright and the 
Hi-bright regions is much more significant than the total quoted uncertainties imply. 
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4. The gas-to-dust ratio 

The spatial resolution of our sub- mm maps and the HI, 12 CO J =1-0 maps from the 
literature permits us to estimate the gas/dust ratio as a function of position in this galaxy. 
Combining the standard expressions for the HI, H 2 and dust mass the gas/dust ratio in 
astrophysical units is expressed as follows: 



— £ « 1.25 x 10 3 — 

M d \v 



P+3 



^hv/kT 



-1 



Shi + 10 X co Sco 



(2) 



where Shi and Sco are the velocity-integrated flux densities of the HI hyperfme transition 
and the 12 CO J=l-0 line respectively in Jy km s _1 , X co is the standard galactic 
conversion factor in units of M (K km s~ x pc 2 ) -1 and S u is the sub-mm flux density at 
frequency v in Jy. We adopted the emissivity law given by Hildebrand (1983), namely 
k(v) = k(i/ D ) {y/v Y, where k(z/ ) = 10 cm~ 2 gr and u = 1196 GHz (250 /zm). 

For a galactocentric radius of r = 20", where most of the FIR emission arises (Telesco 
et al. 1984), we estimate S co = (2800 ± 800) Jy km s" 1 (from the 12 CO J=l-0 channel 
maps, Heifer & Blitz 1995), S ffl = (0.82 ± 0.20) Jy km s^ 1 and S 450 = (8.74 ± 1.95) Jy, 
S 850 = (1.23 ± 0.13) Jy. For (3 = 2 and X co = 5 M (K km s" 1 pc 2 )" 1 (i.e. the standard 
galactic value), the observed 850 /im flux density implies 



Mg 
M d 



315 ( e i6.8/T _ j 



-i 



(3) 



For T = 25 K yields M g /Md ~ 330, similar to the molecular gas-to-dust ratio found in 
spirals using IRAS 100 /im and 60 /im fluxes ( | Young et al. 1986| ; |Young et al. 1989] ). 



Assuming that the actual gas-to-dust ratio in NGC 1068 is similar to that for the Milky 
Way, namely 100-150, we conclude that the application of equation (3) leads to a significant 
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overestimate of this ratio for the inner region (r < 20") of NGC 1068. There is indeed 
independent evidence (e.g. Downes & Solomon 1998) which suggests that the standard 
galactic conversion factor applied to starburst nuclei leads to significant overestimates 
the H 2 content in these regions. In addition, Papadopoulos & Seaquist (1998) provide 
indications that M(H 2 ) in the inner region of NGC 1068 is overestimated by a factor of two, 
which would then produce agreement between the ratio for the nuclear region NGC 1068 
and the Milky way. 

For r > 20", no bright CO emission is detected (see Fig. 2, Heifer & Blitz 1995). 
There is, however, evidence for faint CO emission. By integrating over the area 
20" < r < 60" we find Sco = (1570 ± 470) Jy km s _1 . Over the same area, we estimate 
S(HI) = (16 ± 1) Jy km s" 1 and S 450 = (4.0 ± 0.9) Jy, S 850 = (0.85 ± 0.09) Jy. Using these 
figures, together with equation (2) and the standard value for X co , we obtain, 



For T ~ 10 - 15 K we obtain M g /M d = 70 - 150 and M(H 2 )/M(HI) ~ 5. The good 
agreement of the gas-to-dust ratio with that of the Milky way suggests that most of the 
gas/dust mass is accounted for. 

However, the CO becomes exceedingly faint and essentially undetectable at r > 30". 
If we assume the lowest gas/dust temperature of T ~ 10 K allowed by the S 4 5o/S 8 5o ratio, 
the inferred gas-to-dust ratio for 30" < r < 60" would be ~ 40 with M(H 2 )/M(HI) ~ 1 
(using the standard X co ). Thus, the use of the standard conversion factor in this region 
may underestimate M(H 2 ) by a factor of ~ 5. It seems possible that this molecular gas 
phase is not CO-bright mainly because of low densities (n(H 2 ) < 10 3 cm~ 3 ) rather than low 
gas/dust temperatures. A thermalized 12 CO J=l— line would remain luminous as long 
as T > AEio/k ~5K, and lower temperatures would yield S 4 5o/S 8 5o < 113 which is not 




(4) 
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observed. Assuming similar spatial and velocity filling factors, we estimate that the average 
12 C0 J=l-0 brightness temperature for the outer regions T b ~ 0.07 x T b (r < 20"). For the 
various molecular clouds (size ~ 200 pc) within the r ~ 20" radius T b (r < 20") ~ 10 — 20 K 
(e.g. Planesas, Scoville & Myers 1991), hence for similar size molecular clouds in the 
Hi-bright regions it would be T b ~ 1 K. This is significantly lower than the gas/dust 
temperature inferred by the S 4 5o/S 8 5o ratio over the same regions, thus implying 
sub-thermally excited 12 CO J=l-0. 

These characteristics of the diffuse H 2 phase makes sub-mm measurements very 
valuable since, while CO is very faint and may underestimate the H 2 mass, the dust/gas 
temperature is still high enough to allow sub-mm imaging to reveal the distribution and 
mass of the H 2 mixed with HI, under the assumption of a canonical gas-to-dust ratio. 
We expect this gas phase to be a general feature in spiral galaxies and sensitive sub-mm 
imaging of the regions with high HI column densities is a new tool to study its properties 
through the associated dust emission. 



5. Conclusions 

The results of deep sub-mm imaging of the archetypal Seyfert 2/starburst galaxy 
NGC 1068 can be summarized as follows: 

1. Extended sub-mm emission due to dust with a temperature ofT~10 — 15Kis 
found associated with the regions with the highest HI column densities. A more spatially 
concentrated and warmer (T ~ 20 — 30 K) component is found in the inner starburst 
region where most of the gas is molecular. The estimated gas-to-dust ratio for the inner 
region is M g /Md ~ 330 while for the dust emission associated with the HI gas we find 
M g /M d ~ 70 — 150 depending on the value of the gas/dust temperature. 
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2. Under the assumption that a Milky Way value of M g /M d ~ 100 — 150 applies also 
to NGC 1068, we conclude that the high value of M g /M d in its starburst region is a result 
of an overestimate of the molecular gas mass. This seems to be a systematic effect of using 
the standard galactic conversion factor to convert 12 CO J=l-0 luminosity to molecular gas 
mass in starburst environments. 

3. The diffuse H 2 gas that is mixed with the HI gas is associated with sub-thermally 
excited and very faint CO emission. Furthermore it is more massive than HI, with 
M(H 2 )/M(HI) ~ 5 and the standard galactic conversion factor may underestimate M(H 2 ) 
at low gas/dust temperatures of T ~ 10 K. For this gas phase sensitive sub-mm imaging 
provides a better tool than CO imaging in revealing its distribution and mass. 

We thank Rob Ivison and Lome Avery for conducting a wonderful set of observations 
with SCUBA, and Wayne Holland for friendly and useful advice. We are grateful to Elias 
Brinks for the HI map and Tamara Heifer for the 12 CO J=l map. E. R. S acknowledges 
the support of a research grant from the Natural Sciences and Engineering Research Council 
of Canada. 
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Fig. 1. — Grey scale: Velocity-integrated HI brightness with a range of 100 — 830 mJy 
beam -1 km s -1 , and <7 rms ~ 100 mJy beam -1 km s _1 . 

Contours : 850 fim brightness (top), 450 jum brightness (bottom). The contours are: 

(1, 2, 3, 4, 5, 7, 10, 15, 20, 25, 30, 35)xcr, where cr rms (850 fim) = 10 mJy beam -1 and 

<r rms (450 fim) = 90 mJy beam -1 . The resolution of all the maps is ~ 15". 

Fig. 2. — Gray scale: Velocity-integrated 12 CO J=l-0 brightness from Heifer & Blitz 1995, 
convolved to a resolution of ~ 9", and a range of 100 — 307 Jy beam -1 km s -1 . 
Contours : The bright 450 /im emission at a resolution of ~ 9", the contours are: (8, 10, 12, 
14, 16, 18, 20, 22, 24, 26, 28, 30) xa, where a = 50 mJy beam -1 . 
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Table 1 
Calibration data 



Wavelength #hpbw 


S to t(cal) a 


G, 5G/G b 


a c 

" rms 


(^20, f6o) d 


850 fim 15.25" 


4.56 ±0.17 


265, 15% 


10 


(1.03, 1.30) 


450 /iin 8.75" 


11.5 ± 1.2 


670, 30% 


50 


(1.45, 2.45) 



a The total flux of CRL 618 in Jy, taken from the JCMT secondary calibrators list, 
Sandell 1998. 

b The flux gain in Jy/beam Volts" 1 and its fractional uncertainty estimated from an 
extensive set of beam maps and photometry on CRL 618 and Mars. 

The thermal rms error of the final maps in m Jy/beam. 

d The flux correction factor for circular area with r = 20" and r = 60" radius, estimated 
from beam maps of CRL 618 (see Sandell 1998). If S is the integrated flux density 
within that area, the error-beam corrected flux is S c = S/f. 




Figure 2 



